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In a previous work we have shown that a general two-Higgs-doublet model (THDM) with a very 
light CP-odd scalar can be compatible with electroweak precision data, such as the p parameter, 
BR(& — > sj), Rb, Ab, BR(T — > A7), BR(r; —* A-y), and (g — 2) of muon. Prompted by the recent 
significant change in the theoretical status of the latter observable, we comment on the consequences 
for this model and update the allowed parameter region. It is found that the presence of a very light 
(N ■ scalar with a mass of 0.2 GeV is still compatible with the new theoretical prediction of the muon 

' anomalous magnetic moment. 
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The possibility of a Higgs boson decaying into a pair of light CP-odd scalars was considered in Ref. jj] . Although 
it is very unlikely that this particle can be accommodated in the minimal supersymmetric standard model (MSSM), 
in the light of the restrictions imposed by the current low-energy data on the parameters of this model, a very light 
CP-odd scalar A can still arise in some other extensions of the standard model (SM), such as the minimal composite 
' Higgs model Q, or the next-to-minimal supersymmetric model ||. Therefore, the existence of a very light CP- 
odd scalar not only proves new physics but also casts the most commonly studied MSSM in doubt. Furthermore, 
studying the couplings of the light CP-odd scalar to the SM fermions may help discriminating models of electroweak 
' symmetry breaking - either a weakly interacting model (e.g., the next-to-minimal supersymmetric model) or a strongly 
interacting model (e.g., the minimal composite Higgs model). Apart from the above implications arising from the 
' existence of a very light CP-odd scalar, our main interest in studying this particle stems from the fact that its 
phenomenology is indeed rather exciting: an interesting aspect of a light A is that if its mass Ma is less than twice 
O ' that of the muon m^, i.e. less than about 0.2 GeV, it can only decay into a pair of electrons (A — > e + e~) or photons 
(A — ► 77). Hence, the decay branching ratio BR(A — > 77) can be sizable. Consequently, A can behave like a 
fermiophobic CP-odd scalar and predominantly decay into a photon pair, which would register in detectors of high 
;> energy collider experiments as a single photon signature when the momentum of A is much larger than its mass jlj . 
T ■ In a previous work |ij we performed an extensive analysis within the framework of the two Higgs doublet model 
(THDM) and found that a very light CP-odd scalar can still be compatible with precision data, such as the p 
parameter, BR(6 — > S7), Rb, Ab, BR(T — » A7), and the muon anomalous magnetic moment a M j|. We considered 
different values for sin 2 (/? — a) and found the constraints imposed on the remaining parameters of the model, which 
we summarize in Table Q, where Mg (M/,) stands for the heavy (light) CP-even scalar and M H ± for the charged 
scalar. As for the soft breaking term /112, it is not involved in any of the above observables, so they cannot be used to 
constrain it. Since has no relevance for the present discussion (the purpose of this work is to update the bounds 
derived from the changes in the status of the theoretical value of the muon anomaly) , we refer the reader to Ref. [Q , 
where LEP-2 data were used to set bounds on this parameter. In obtaining the bounds shown in Table | we have 
used the lower value of 110 GeV for M%. We recall that the LEP-2 direct search bound requires Mh > 114.1 GeV at 
the 95% GL. ||. However, in the presence of new physics such a bound can be substantially relaxed. As explained 
in Ref. the reason why the LEP-2 bound (Mh > 114.1 GeV) does not apply in our model is because this bound 
is based on the SM specific value of BR(/i — ► 66) ||. In the THDM, the new h — > AA decay mode can significantly 
reduce the h — > bb branching ratio. This was clearly illustrated in the Fig. 9 of Ref. B for some allowed parameter 
space of the model. In any case, the new decay channel (h — > AA) registers as a di-photon signature (h — > 77) for 
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TABLE I: Constraints from the low energy data for Type-I and Type-II THDMs, with M A = 0.2 GeV. The old calculation of 
the hadronic light by light contribution to a M , cf. Table |n|, was used together with the calculation of Ref. [^] for the hadronic 
vacuum polarization. When sin 2 ((3 — a) = 1 there is no Mh dependence on p, otherwise, we assume Mh = 110 GeV. 

Constraint Type-I THDM Type-II THDM 

(g-2)f, tan/3>0.4 tan /3 < 2.6 

[tan $ > 1] b -> S7 M H + > 100 GeV M H+ > 200 GeV 
[0.5 < tan/3 < 1] b -> s-y M H+ > 200 - 350 GeV 
[0.6 < tan/3 < 1] R b M H+ > 200 - 600 GeV M H + > 200 - 600 GeV 
[sin 2 (/3 - a) = 1] Ap M H ~ M H + M H ~ M H + 

[sin 2 (/3 - a) = 0.8] Ap M H ~ 1.2 M ff+ Mh ~ 1.2 M H + 

[sin 2 (/3 - a) = 0.5] Ap M H ~ 1.7 M H+ M g ~ 1.7 M H + 



which LEP-2 has already set a lower bound. By taking both the AA and bb decay modes into consideration, a lower 
bound for Mh > 103 GeV can be established in our light A scenario 0]. 

At this point we would like to emphasize that, given the recent measurements of at Brookhaven National 
Laboratory (BNL) M, the bounds on new physics effects imposed by the muon (g — 2) data depend largely on the 
theoretical value predicted by the SM for the nonperturbative hadronic contribution to a^. In our analysis 01 , we 
followed a conservative approach and considered various_ predictions for the hadronic correction a„ ad || ^, |lC| [ll] , 



which in fact has been the source of debate recently [[L2|, 13, |LJ]. For instance, the bounds shown in Table | were 
obtained from the calculation presented in Ref. || , which was the one allowing the largest parameter space. 

After the completion of our work, it was evident that the latest precision measurement of a M at BNL JtJ along 
with some theoretical predictions for ajj ad disfavored the presence of a light A in the THDM. As is well known, the 
BNL data opened the prospect for new physics as the experimental value of a M appeared to be more than 2.6 a 
above the theory prediction based on some calculations of the hadronic vacuum polarization. At the one-loop order, 
a light CP-odd scalar can give a significant negative contribution to a M , making it harder for this type of model to 
be consistent with experiment. However, the two-loop calculation can yield a large correction to the one-loop result 
as pointed out in jlq ]. Although this fact seems to contradict perturbation theory, the unusual situation in which a 
two-loop diagram can give a contribution of similar size or even larger than that from the one-loop diagrams within 
a perturbative calculation was noted first by Bjorken and Weinberg when evaluating the Higgs scalar contribution 
to the /j, 



e7 decay |16 



It is straightforward to see that this situation also occurs in the calculation of the Higgs 
scalar contribution to a^. The reason is that the coupling of the Higgs scalar to the muon enters twice in the one- loop 
diagram, whereas at the two-loop level there appears a diagram in which this coupling enters just once, together 
with a line where the Higgs scalar couples to a heavy fermion pair (see Fig. 0). This gives rise to an enhancement 
factor, due to the couplings, that compensates the suppression factor g 2 /(167r 2 ), due to an additional loop. It turns 
out that the diagram of Fig. [j] (c) gives by far the most dominant contribution at the two-loop level. Therefore, we 
do not expect large uncertainties arising from unknown higher order terms. In our previous analysis, even when we 
considered the two-loop calculation for the CP-odd scalar contribution to a M , together with the hadronic correction 
quoted in Ref. Q, i.e. that by Davier and Hocker j poj , we found that there was no allowed parameter space (in the 
type-II THDM) in the tan/3 vs. Ma plane when Ma was below 3 GeV. Nevertheless, there were other SM calculations 
yielding close enough to the experimental value as to allow a very light A. 
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FIG. 1: Contribution from the THDM to the anomalous magnetic moment of muon: (a) neutral Higgs bosons, (b) charged 
Higgs boson, and (c) the leading two-loop contribution from the CP-odd scalar. 
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TABLE II: Contributions to the anomalous magnetic moment of muon in the SM JlSj], prior to the discovery of a wrong sign in 
the pion pole correction to ajj ad (l.b.l), which significantly changed the a^ hoory prediction. All values are given in units of 10 -11 . 

Contribution SM prediction 



ajf 116584705.7 (1.8) 

a£ ad (l.b.l) -79.2 (15.4) ° 



.. h ad 



(h.o.) -101 (6) 



"This value has been found to be wrong in Ref. Jl?} ]. 



Since the publication of |7]] , there has been a lot of controversy regarding the theoretical value of the muon anomalous 
magnetic moment. It is evident that before claiming the presence of any new physics effect, an extensive reexamination 
of every contribution to a M is necessary [[l2| [u^, |l4| . Along these lines, a reevaluation |T^] of the hadronic light by light 
contribution to a M found a sign error in earlier calculations |l8| ] of this contribution, which has resulted in a significant 
change of the a M prediction. Once the corrected value is taken into account, the discrepancy between experiment and 
theory reduces down to the level of 1.6 a. In the light of this result, we believe it is worth revisiting our work and 
reexamining our previous bounds. 



II. ALLOWED PARAMETER RANGE FOR M A AND tan/3 

The SM prediction of a M is composed of the following three parts |19| : 



a theory = a QED + a woak + ^ad (!) 

where the electroweak corrections have been computed to a very good accuracy: they have a combined error of 
the order of 5 x 10 -11 , which is already about one order of magnitude smaller than the ultimate goal of the E821 
experiment [Q. In contrast, the hadronic contribution a^ ad contains the bulk of the theoretical error (~ 70 x 10~ n ) 
and can be decomposed into three parts [fl9|| , namely the hadronic vacuum polarization contribution a M (h.v.p.), the 
hadronic light-by-light correction a M (l.b.l), and other hadronic higher order terms a Al (h.o.): 



a* ad = a^h.v.p.) + a;r d (l.b.l) + aj^h.o.). (2) 

In our previous analysis we used the values shown in Table [II] for each contribution to aJ; hcory |nj together with the 
aj 1 1 ad (h.v.p.) predictions to be discussed below. In the months following the publication of our work, a new situation 
arose: the sign of the pion pole contribution to the hadronic light by light correction was found to be wrong |l7fl . Very 
interestingly, this contribution alone represents about 70% of the full a^ ad (l.b.l.). It turns out that after correcting this 
mistake, the ajj ad (l.b.l.) value gets significantly changed and even its sign gets flipped. As a result, the discrepancy 
between the experiment and theory reduces down to the level of 1.6 a. Subsequent publications have confirmed this 

we list the most recent evaluations of a^ ad (l.b.L). In addition, there is one more 
calculation that is based on chiral perturbation theory |23[: 



finding POL ElL B2l . In Tabic [II 



a^ ad (l.b.l.) = ( 55±™ + 31C ) 



x 10" 



where C is an unknown low-energy constant that parametrizes some subdominant terms. We will not consider this 
result here but only mention it as an example of a calculation that is still open to debate. For the purpose of this 



work we will take an average of the top three results shown in Table [II and study the consequences on the allowed 
parameter space of the THDM. 

After introducing the corrected value of et^ ad (l.b.l.), the sum of all the contributions to a^ hcory except a{j ad (h.v.p.) 

is: 



a thcory _ a had( h v p ) = 11658 4 8 45.3 ( 17 1 ) x 1( )-U 



(3) 
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TABLE III: The most recent evaluations of the hadronic light by light contribution to (l.b.1. 



contrast with the wrong one shown in Table [fj|. 



These corrected values 



Author 


aj; ad (l.b.l.) x 10 11 


KN 111 


83(12) 


hkH 


89(15) 


BPPMI 


83 (32) 


BCM p>| 


56 a 



"This value accounts only for the pion pole contribution. 



TABLE IV: Some of the most recent calculations of ajl ad (h.v.p.) together with the respective theory prediction a* hoory and the 
discrepancy Aa M between experiment and theory. The last column represents the bounds on any new physics contribution a^ p 
at the 95% C.L. All of the values are given in units of 10 -11 . 



Author 


ar>-v.p.) 


theory a 


Allowed range for 


adhJI 


7011(94) 


116591856.3 (95.54) 166.7(179.53) 


[-185 


-519] 


DH 




6924(62) 


116591769.3 (64.31) 253.7(164.92) 


[-70- 


577] 


J 


iff 


6974(105) 


116591833.3 (112.3) 189.7(188.8) 


[-180 


-560] 


N 




7031(77) 


116591876.3 (78.88) 146.7(171.25) 


[-189 


-482] 


TY 




6952(64) 


116591797.3 (66.25) 225.7(165.81) 


[-99- 


551] 



where the errors have been composed quadratically. 1 As for the et M (h.v.p) term, its evaluation has also been the 
source of renewed interest lately. 2 As in our previous work, here we will use a conservative approach and consider 
some representative evaluations of a M (h.v.p). In the second column of Table |y| we show some of the most recent 
results, which were compiled in [p4[ , whereas in the third column we show the full theory prediction, which is obtained 
after adding up each value in the second column to Eq. (||). 

As for the experimental value a° xp , the data obtained during the 1999 running period combined with previous 
measurements give M 



a c * p = 116592023 (152) x 10" 



(4) 



One thus can obtain the discrepancy between experiment and theory Aa M 



- a' hoory for each different evaluation 



of a A1 (h.v.p), as shown in the fourth column of Table [V. Finally, if we assume that the discrepancy between theory 
and experiment is to be ascribed to new physics effects, we can obtain the bounds shown in the last column of the 
same Table for the new physics contribution to the anomalous magnetic moment at the 95 % C.L., which is denoted 
by a^ p . Those bounds on a^ p should be compared to those used in our previous analysis, cf. Eq. (2) in Ref. Qj. 

Given the new bounds on a^ p , we update the constraint imposed by it on the tan (3 — Ma plane within the THDM. 
The analytical expressions for the contribution of either a CP-even or a CP-odd scalar (Fig. |) can be found in 
Appendix A of Ref. Q. We will use the two-loop calculation for the contribution from the CP-odd scalar [ p!5[ . In 
order to satisfy the bounds shown in Table |, we are assuming that the remaining four Higgs scalars are much heavier 
than the CP-odd scalar A, so their contribution to turns out to be negligibly small as compared to that coming 
from the latter. Also, we are considering that sin 2 (/3 — a) — 1. The reason why we make this choice is because in our 
scenario with a very light CP-odd scalar the most convenient way to meet the constraint imposed by the p parameter 



is to have Mh and M H + nearly degenerate and sin 2 (/3 — a) close to 1 
the bounds arisen from the theoretical predictions based on the DH fllO 



-.had 



(h.v.p 



For comparison purposes, we will analyze 
J PJ, N |2~1 and TY || calculations of 



which are the most representative and recent ones. We would like to note that, as observed through Fig. 



to Fig. |6, the bounds from the J [|11| and N [|25| calculations are almost indistinguishable. 

In Figs. ^-[| we show the allowed regions in the tan/3 — Ma plane for both types of TEIDMs. In Fig. |2|, which shows 
the low Ma regime, it can be seen clearly that even if one considers the DH calculation of a^ ad , which is the one with 
the smallest error, there is still the possibility of having a CP-odd scalar with a mass of the order of 0.2 GeV in the 
type-II THDM as long as tan/3 < 1.43, whereas for a type-I THDM tan/3 has to be greater than 0.87. This is a very 
significant change with respect to the results obtained when using the old (uncorrected) value of a^ hcory . In that case, 



1 Throughout this work we will systematically compose the errors in quadrature. 

2 For a summary of the most recent evaluations of ajj ad , see Refs. [ p4{ and pi) ]. 
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the DH calculation did not allow for a light CP-odd scalar in either type of THDM, though other calculations did 
allow such a possibility. 




FIG. 2: The regions (above the curves for type-I and below the curves for type-II THDM) in the tan/3 versus Ma plane allowed 
by the a M data at the 95% CL. Four different curves are displayed depending whether the SM prediction is obtained from the 
DH, J, N or TY calculation of o)i ad (h.v.p.). The two-loop contribution from the light A has been used. 



□H ■ 
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TY ■ 



0.1 ' ' — 1 — 1 — 1 — 

0.01 0.1 1 10 100 

M A [GeV] 

FIG. 3: The region (above the curves) in the tan/3 versus Ma plane of a type-I THDM allowed by the a M data at the 95% 
CL. The allowed regions based on the DH, J, N and TY calculations are above the curves. The two-loop contribution from the 
light A has been used. 

As stated above, so far our results have been derived from the two-loop contribution from the CP-odd scalar to aj^ p . 
It is also interesting to repeat the above analysis using only the one-loop calculation for a^ p . Its result is depicted 
in Figs. ^ and ^. The old (uncorrected) theory prediction based on the DH calculation required any new physics 
contribution to a M to be positive. However, the one-loop contribution from a light CP-odd scalar is always negative. 
Therefore, the old SM theory prediction for combined with the THDM one-loop correction strongly disfavored the 
existence of a very light CP-odd scalar. This is to be contrasted with the conclusion drawn from the corrected value 
of a^ llcory . In that case, there is indeed an allowed region of tan/3 when Ma ~ 0.2 GeV, though this region is smaller 
than the one allowed by the two-loop calculation of a^ p (cf. Figs. 3 and 5, and Figs. 4 and 6). As shown in Fig. |J, 
there is an interesting feature in the tan/3 versus Ma plane of a type-II THDM when Ma is around 2.6 GeV. It is 
because for Ma ~ 2.6 GeV, the two- loop contribution from a light CP-odd scalar becomes as large as the respective 
one-loop contribution but with an opposite sign, so the total effect cancels. 



A. Bounds on tan/3 from meson decays 



For completeness we now turn to analyze the bounds obtained on THDMs with a very light CP-odd scalar from 
meson decays. A very light Higgs scalar (CP-odd or CP-even) can be a decay product of some hadrons, like the ij 
and T mesons. For the latter, a measured upper bound to the X + 7 decay channel has been set p^j that can be 
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M A [GeV] 

FIG. 4: The regions (below the curves) in the tan/3 versus Ma plane of a type-II THDM allowed by the a M data at the 95% 
CL. The allowed regions based on the DH, J, N and TY calculations are below the curves. The two-loop contribution from the 
light A has been used. 




FIG. 5: Same as Fig. H, but only the one-loop contribution from the light A is considered. 




M A [GeV] 

FIG. 6: Same as Fig. rt, but only the one-loop contribution from the light A is considered. 



used to constrain the Abb coupling. Denote the Yukawa coupling of A b b to be kdmb/v, with kd = tan/3 (cot /3) in 
the type-II (type-I) model. Then, the data of the meson decay T — > 7 + X requires kd < 1. (We refer the reader to 
Refs. [[| [| for a detailed discussion.) 

As shown in Ref. |2§[| , there is another decay process that can strongly constrain tan/3, namely 77 — ► irS, where S 
is a very light CP-even scalar. Those results can be translated into the case of a CP-odd scalar. In particular, the 
experimental upper limit 
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BR(ry -» 7r°e + e~) < 5 x 10~ 5 (5) 

can be used to obtain the following constraint on a THDM CP-odd scalar with mass Ma lying in the range 2m e < 
M A < 2m M : 

(k d -k u ) 2 \(l ,4 .^r) <l-5 (6) 

V ™ 2 ™ 2 y 

where fe u is cot/3 for either type-I or type-II THDM and k d has been defined above. The function A is given by 
A 2 (a, b, c) = a 2 + b 2 + c? — 2ab — 2ac — 2bc. From here we can conclude that cot/3 > 0.65 for type-I THDM and 
0.55 < tan/3 < 1.8 for type-II THDM. We thus can confirm that the hadron decay data together with the muon 
(<7 — 2) measurement require tan /3 to be of order 1 if there exists a very light pseudoscalar with a mass smaller than 
2m M . 



III. OVERALL DESCRIPTION OF THE GENERAL THDM WITH A LIGHT A 



Once the allowed parameter range for tan (3 and Ma has been updated, there remains five other parameters to 
consider: the CP-even neutral Higgs mixing angle a, the soft breaking term /ii2 and the three other Higgs masses: 
M/j, Mh and M H +. Since we already know that tan/3 has to be of order 1 we can address the status of the charged 
Higgs mass M H + independently of the other parameters. It turns out that both the b — > S7 and the Rb data require 
H + to be considerably heavy 0, 29 : 



M H + > 350 GeV . (7) 

Such a high lower bound for the H + mass affects the allowed values of the mixing angle a. In Ref. Q we show 
that the p parameter requires Mh and M H + to be very correlated depending on the value of sin 2 (/3 — a). In fact, 
if a very light CP-odd scalar is to be allowed, the easiest way to satisfy the bound imposed by p ~ 1 is to have 
Mh and M H + degenerate and sin 2 (/3 — a) = 1. With this choice, Mh is not restricted since it does not contribute 
to the p parameter. As we consider values of sin 2 (/3 — a) smaller than 1, it turns out that p is very sensitive to the 
masses of H and H + . For instance, if sin 2 (/3 — a) = 0.5, Mh must be at least of the order of 500 GeV Qj. Generally 
speaking, our conclusion on the bounds on a very light CP-odd scalar in the THDM does not change significantly for 
0.5 < sin 2 (/3 — a) < 1 as long as the other Higgs bosons in the model are heavy enough. For a very small value of 
sin 2 (/3 — a), much less than 0.5, the p-parameter data would have required the mass of H to be at the TeV order. 



IV. CONCLUSION 



In conclusion, with the recent correction to the SM prediction of a^, the current muon (g — 2) data, together with 
other precision data (cf. Table I), still allows a light (Ma ~ 0.2 GeV) CP-odd scalar boson in the THDM. Due to this 
new development in the SM theory calculation of muon's (g — 2) , the allowed range of tan /3 in the Type-I or Type-II 
THDM is modified, and our result is summarized in Table |v| It is interesting to note that the phenomenology at 
high energy colliders predicted by the THDM with a light CP-odd Higgs boson is dramatically different from that 
predicted by the usual THDM in which the mass of the CP-odd scalar is at the weak scale. A detailed discussion 
on this point can be found in Ref. Q. In particular, various potential discovery modes were studied in there: it was 
found that the Fermilab Tevatron, the CERN large hadron collider (LHC) and the planned e + e~ linear collider (LC) 
have a great potential to either detect or exclude a very light A in the THDM. 

Finally, we note that while a light CP-odd scalar in THDM is still compatible with all the precision data, it has been 
shown recently in Ref. |3(| that a light CP-odd scalar in the MSSM will violate the constraint derived from the Zbb 
coupling. This is because in the MSSM, the masses of the five Higgs bosons are related by the mass relations required 
by supersymmetry. Hence, with a light CP-odd scalar, the mass of the other Higgs bosons cannot be arbitrary large, 
and it is difficult to yield the decoupling limit when calculating low energy observables. 



Note added: 
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TABLE V: Constraints on tan f3 from the muon (g - 2) data for Type-I and Type-II THDM, with M A = 0.2 GeV, based on 
various SM theory predictions of a^ ad (h.v.p). The two-loop contribution for the CP-odd scalar has been used. 

Theory prediction Type-I THDM Type-II THDM 





tan/3 > 0.87 


tan/3 < 1.43 




tan/3 > 0.54 


tan/3 < 2.19 




tan/3 > 0.53 


tan/3 < 2.24 




tan/3 > 0.73 


tan/3 < 1.67 



During the review process of this manuscript, the muon (<? — 2) collaboration announced a new result based on 
data collected in the year 2000 in which the experimental uncertainty has been reduced to one half that of 
the previous measurement while the central value of a® xp remains about the same. [The new data yields a® xp = 
11659204(7) (5) x 10~ 10 (0.7 ppm).] According to the latest experimental data, we have updated Figs. || to ^| in 
this paper to Figs. to y. The new data suggests that a very light CP-odd scalar is not allowed in the type-I or 
type-II THDM based on the SM calculation done by DH and TY §|. However, based on the N || and J 
calculations, a very light CP-odd scalar is still possible though the allowed parameter space of the THDM has been 
tightly constrained. 
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FIG. 7: Same as Fig. ^, but with the latest experimental data from the muon (g — 2) collaboration 
region in this range of parameters according to the DH fllCI and TY [M calculations. 



There is no allowed 



M A [GeV] 



FIG. 8: Same as Fig. ^, but with the latest experimental data from the muon (g — 2) collaboration pjj . There is no allowed 
region in this range of parameters according to the DH |1CJ and TY |E4| calculations. 
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TY 
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M A [GeV] 



FIG. 9: Same as Fig. but with the latest experimental data from the muon (g — 2) collaboration | p]| . The region allowed 
by the DH Q and TY J^|] calculations is bounded by the respective lines. 



M fl [QeV] 



FIG. 10: Same as Fig. ^, but with the latest experimental data from the muon (g — 2) collaboration pjj . There is no allowed 
region in this range of parameters according to the DH ||10| and TY |M calculations. 



M fl [GeV] 



FIG. 11: Same as Fig. |6| but with the latest experimental data from the muon {a — 2) collaboration considered There is 
no allowed region in this range of parameters according to the DH |10| and TY |24[ calculations. 
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